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Abstract. Motivated by recent analytical and numerical work on two- and 
three-dimensional convection with imposed spatial periodicity, we analyse three 
examples of bifurcations from a continuous group orbit of spatio-temporally 
symmetric periodic solutions of partial differential equations. Our approach 
is based on centre manifold reduction for maps, and is in the spirit of earlier 
work by Iooss (1986) on bifurcations of group orbits of spatially symmetric 
equilibria. Two examples, two-dimensional pulsating waves (PW) and three- 
dimensional alternating pulsating waves (APW), have discrete spatio-temporal 
symmetries characterized by the cyclic groups Z n , n = 2 (PW) and n = 4 
(APW). These symmetries force the Poincare return map A4 to be the n th 
iterate of a map Q: A4 = Q n . The group orbits of PW and APW are generated 
by translations in the horizontal directions and correspond to a circle and a 
two-torus, respectively. An instability of pulsating waves can lead to solutions 
that drift along the group orbit, while bifurcations with Floquet multiplier +1 of 
alternating pulsating waves do not lead to drifting solutions. The third example 
we consider, alternating rolls, has the spatio-temporal symmetry of alternating 
pulsating waves as well as being invariant under reflections in two vertical planes. 
When the bifurcation breaks these reflections, the map Q has a "two-symmetry," 
as analysed by Lamb (1996). This leads to a doubling of the marginal Floquet 
multiplier and the possibility of bifurcation to two distinct types of drifting 
solutions. 
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1. Introduction 

Techniques for analysing symmetry-breaking bifurcations of T-invariant equilibria of T- 
equivariant differential equations are well-developed in the case of compact Lie groups T 
(Golubitsky et al 1988). The motivation for developing these methods comes, in large part, 
from problems of pattern formation in fluid dynamics (see, for example, Crawford & Knobloch 
1991). In the simplest cases, the symmetry-breaking bifurcation corresponds to a pattern- 
forming instability of a basic state that is both time-independent and fully symmetric, for 
example, a spatially uniform equilibrium solution of the governing equations. A symmetry- 
breaking Hopf bifurcation of this spatially uniform state often leads to time-periodic solutions 
that break the translation invariance of the governing equations and that have spatio-temporal 
and spatial symmetries. In this paper we address bifurcations of such periodic orbits, which 
have broken the translation invariance but have retained a discrete group of spatio-temporal 
symmetries. 

We consider problems posed with periodic boundary conditions, for which there is an 
S 1 symmetry associated with each direction of imposed periodicity. If this symmetry is 
broken by an equilibrium solution, then the solution is not isolated; there is a continuous 
family of equilibria related through the translations. An instability of this solution can excite 
the neutral translations modes(s) and lead to new solutions that drift along the translation 
group orbit. This is the case, for example, in the "parity-breaking bifurcation": a reflection- 
symmetric steady state undergoes a symmetry-breaking bifurcation to a uniformly translating 
solution. Another example of a bifurcation leading to drift has been observed in two- 
dimensional convection: when the vertical mirror plane of symmetry that separates steady 
counter-rotating rolls is broken in a Hopf bifurcation, the resulting solution, called a direction- 
reversing travelling wave or pulsating wave (PW), drifts to and fro (Landsberg & Knobloch 
1991; Matthews et al 1993). This periodic orbit is invariant under the combination of advance 
of half the period in time with a reflection; any drift in one direction in the first half of the 
oscillation is exactly balanced by a drift in the other direction in the second half, so there 
is no net drift during the oscillation. Similarly in thee-dimensional convection with spatial 
periodicity imposed, for example, on a square lattice, a symmetry-breaking Hopf bifurcation 
from steady convection in a square pattern can lead to alternating pulsating waves (APW), 
which are invariant under the combination of advance of one quarter the period and rotation 
by 90° (Rucklidge 1997). These solutions drift alternately along the two horizontal coordinate 
directions, but again have no net drift over the whole period of the oscillation. 

There have been a number of studies of bifurcations of compact group orbits of (relative) 
equilibria. Iooss (1986) developed an approach based on centre manifold reduction to 
investigate bifurcations of Taylor vortices in the Taylor-Couette problem. Specifically, he 
analysed bifurcations in directions orthogonal to the tangent space to the group orbit of 
equilibria, with the neutral translation mode incorporated explicitly in the bifurcation problem. 
Krupa (1990) provided a general setting for investigating bifurcations of relative equilibria 
that focuses on the local dynamics in directions orthogonal to the tangent space to the 
group orbit. He shows that the resulting bifurcation problem is X-equivariant, where £ 
is the isotropy subgroup of symmetries of the relative equilibrium, and, building on work 
of Field (1980), provides a group theoretic method for determining whether or not the 
bifurcating solutions drift. Aston et al (1992), and Amdjadi et al (1997) develop a technique 
for numerically investigating bifurcations of relative equilibria in 0(2)-equivariant partial 
differential equations, and apply their method to the Kuramoto-Sivashinsky equation. Their 
approach isolates one solution on a group orbit, while still keeping track of any constant drift 
along the group orbit. 

In this paper we investigate bifurcations of time-periodic solutions that are not isolated 
as they have broken the translation invariance, but that do possess a discrete group of spatio- 
temporal symmetries. Our approach is similar to that of Iooss (1986). However, we are 
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interested in instabilities of periodic solutions, so we use centre manifold reduction for Poincare 
maps. We are particularly interested in determining whether the symmetries of the basic state 
place any restrictions on the types of bifurcations that occur, and whether the bifurcating 
solutions drift along the underlying group orbit or not. We consider three examples that 
are motivated by numerical studies of convection with periodic boundary conditions in the 
horizontal direction(s). First we investigate bifurcations of the pulsating waves and alternating 
pulsating waves described above. These solutions have discrete spatio-temporal symmetries 
Z 2 and Z 4 , respectively. The group orbit of the pulsating waves is 5 1 , while the group orbit 
of the alternating pulsating waves is a two-torus, due to imposed periodicity in two horizontal 
directions. The third example we treat in this paper is alternating rolls (AR), which have 
the same spatio-temporal symmetry as APW but are also invariant under reflection in two 
orthogonal vertical planes (Silber & Knobloch 1991). 

The Z n (n = 2, 4) spatio-temporal symmetry of the basic state places restrictions on the 
Poincare return map M; specifically, we show that it is the n th iterate of a map Q. A direct 
consequence of this is that period-doubling bifurcations are nongeneric (Swift & Wiesenfeld 
1984). Throughout the paper we restrict our analysis to bifurcation with Floquet multiplier +1; 
we do not consider Hopf bifurcations. We also restrict attention to bifurcations that preserve 
the spatial-periodicity of the basic state. 

Our paper is organized as follows. In the next section we lay the framework for our 
analysis in the setting of a simple example, namely bifurcation of pulsating waves. We show 
how the spatio-temporal symmetry is manifest in the Poincare return map. Section 3 considers 
bifurcation of the three-dimensional analogue of pulsating waves, namely alternating pulsating 
waves. Section 4 considers bifurcations of alternating rolls. For this problem we need to 
consider six different cases, which we classify by the degree to which the spatial, and spatio- 
temporal symmetries are broken. In the case that the spatial reflection symmetries are fully 
broken by the neutral modes, the Floquet multiplier +1 is forced to have multiplicity two, 
and more than one solution branch bifurcates from the basic AR state. In one case we find a 
bifurcation of the AR state leading to two distinct drifting solutions. We present an example of 
one of the drifting patterns that is obtained by numerically integrating the equations of three- 
dimensional compressible magnetoconvection. In the course of the analysis of bifurcations of 
alternating rolls, we make contact with the work on ^-symmetries of Lamb & Quispel (1994) 
and Lamb (1996). Section 5 contains a summary and indicates some directions for future 
work. 

2. Two dimensions: pulsating waves 

We write the partial differential equations (PDEs) for two-dimensional convection symbolically 
as: 



where U represents velocity, temperature, density, etc. as functions of the horizontal 
coordinate x, the vertical coordinate z and time t; [i represents a parameter of the problem; 
and T is a nonlinear operator between suitably chosen function spaces. We assume periodic 
boundary conditions, with spatial period £, in the x-direction. 

The symmetry group of the problem is 0(2), which is the semi-direct product of Z 2 , 
generated by a reflection n x , and an 50(2) group of translations r , which act as 




(1) 



r a : x — ► x + a (mod £), 



(2) 



where r e is the identity and r a The PDEs (1) are equivariant under the action of 

these symmetry operators, so J r (r C/';/x) = T a T(JJ\ii) and T{k x XJ\[i) = k x T{U; //), where r a 
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and k x act on the functions as follows: 

T a U(x,z,t) = U(x-a,z,t), K x U(x,z,t) = M Kx lI(-x,z,t). (3) 

Here M Kx is a matrix representing k x ; it reverses the sign of the horizontal component of 
velocity and leaves all other fields in U unchanged. 

Suppose that when the parameter \i = 0, there is a known pulsating wave 
solution U (x, z, t) of (1) with temporal period T and spatial period A = £/N, where N specifies 
the number of PWs that fit into the periodic box. The symmetries of U are summarized as 
follows: 

U (x,z,t) = K x U (x,z,t+ \T) = U (x,z,t + T) =T X U a (x,z,t). (4) 

There is a continuous group orbit of PWs generated by translations: U a = t u Uq. We are 
interested in bifurcations from this group orbit. Following the approach developed by Iooss 
(1986) and Chossat & Iooss (1994) for studying instabilities of continuous group orbits of 
steady solutions, we expand about the group orbit of periodic solutions as follows: 

U(x, z, t) = T c{t) (U Q (x, z, t) + A(x, z, t)). (5) 

Here translation along the group orbit is given by t c u\ , where c is a coordinate parameterizing 
the group orbit. Small perturbations, orthogonal to the tangent direction of the group orbit, 
are specified by A(x,z,t). The expansion (5) is substituted into the PDEs (1) and, after 
suitable projection that separates translations along the group orbit from the evolution of the 
perturbation orthogonal to it, we obtain equations of the form (see Chossat & Iooss (1994)): 

A A A r 

— = g(A,U ;v), - = h{A,U ;fj), (6) 

where Q and h satisfy Q(0,U o ;0) = and h(0,U ;0) = 0. An important consequence 
of the translation invariance of the original PDEs is that Q and h do not depend on the 
position c along the group orbit; the equation for the drift c is decoupled from the equation 
for the amplitude of the perturbation A. Here we find it convenient to keep track of the 
explicit time dependence of Q and h, which enters through their dependence on the basic 
state Uo, by listing Z7 as one of the arguments of Q and h. We determine how the spatio- 
temporal reflection symmetry of U is manifest in the equations for c and A by noting that if 
T c (t)(Uo(x,z,t) + A(x, z,t)) is a solution of the PDEs (1) , then so is 



Hence 



K x T c{t) (U (x,Z,t) +A(x,z,t)) = T- c{t) (K x U (x,Z,t) + k x A(x, z,t)). (7) 

G(k x A, k x U ;h) = K x Q(A,U ;fj,), 

(8) 

h(n x A, n x U ;n) = -h(A, U ;fi). 

Since our basic state U is T-periodic, we seek a map that gives the perturbation A at 
time t = T given a perturbation A(0) at some initial time t = 0. Specifically, we define 
a time advance map Ml acting on the perturbation A by A(t) = Ml(A(0)). We adopt 
the approach of Swift k. Wiesenfeld (1984) and split the time interval from to T into 
two stages using the symmetry property of the underlying pulsating waves. Specifically, 
since K x A(t) satisfies d ^/^ = G{k x A, k x U ; fi) and k x U q (x, z, t) = U Q (x, z,t + -j), we have 
K x A(t) = MjffK^lO)); hence 

k x MI = M T T % +t k x . (9) 
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Figure 1. Illustration of k x Mo = M^j^* n x - In this example, the reflection k x changes the 
sign of the horizontal coordinate. The PW periodic solution is shown as a dotted line, (a) A 
perturbation at t = is advanced in time by an amount t (the solid line, which stays close 
to the broken line on the periodic orbit), then the system is reflected, (b) We arrive at the 
same final position if we reflect (so now the perturbation is about the PW at t = |T) and then 
advance in time by the same amount. 



Advancing the perturbation by a time t starting from time and then reflecting the whole 
system is equivalent to reflecting the whole system then advancing by a time t starting from 
time \T (see Figure 1). It follows immediately that the full period map Mq can be written 
as the second iterate of a map Q: 



Ml 



Ml /2 nlMl /2 = (n x Ml /2 ) 2 = 



(10) 



Rather than consider the full period map Mq, we will consider the map Q = n x M^ 2 . The 
map Q has no special property under reflections, but it commutes with translations t\, which 
leave the underlying pulsating waves invariant: Qt\ = T\Q. 

The dynamics of the perturbation is now given by the map Q: A n+ i = Q(A n ;/J.), where 
each iterate corresponds to advancing in time by \T and reflecting; thus A(^T) = k x Ai, 
starting from Aq at time 0. In order to compute the drift C\ of the solution at time |T, 
we integrate the dc/dt equation (6) for a time |T, starting at a position c and with initial 
perturbation A(0) = A : 



f T / 2 

ci = c + / h(Ml(A ), U (t); fi) dt = c + h(A ; n). 
Jo 

Then, after a second half-period, 

c 2 = Cl +/ fc(M^ /2 (A(!T)),L/o(*);M)dt 

JT/2 

= ci+ f h(MT /2 {K x Ai),U (ty,n)dt 

JT/2 

= ci + / hi^M^^iA^, K x U {t - T/2);n) dt 



(11) 



T/2 
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,T/2 

= Cl -/ h(M\{Ai),U (i!y,ii)ti! 
Jo 

= d-h(Ai;n). (12) 

Thus the combined dynamics of the perturbation and translation can be written as 

A n+1 = G(A n ;n), c n +i = c n + (-l) n h(An-,n). (13) 

Since the unperturbed PW is a non-drifting solution of the problem at fi = we have 
£7(0; 0) = and h(0;0) = 0. Moreover, the spatial periodicity of U places some symmetry 
restrictions on Q and h; specifically, Q(t\A; fi) = T\Q{A; fi) and Ii(t\A; fi) = h(A; fi). 

We turn now to the codimension-one bifurcations of the PW, which are the trivial fixed 
points A = 0, c = c of (13) when fi = 0. The map (13) always has one Floquet multiplier 
(FM) equal to one because of the translation invariance of the c part of the map. Bifurcations 
occur when a FM of the linearization of Q crosses the unit circle: either a FM = 1, or a 
FM = —1, or there is a pair of complex conjugate FMs with unit modulus. Because we have 
assumed periodic boundary conditions in the original PDEs, we expect the spectrum of the 
linearization to be discrete and the centre manifold theorem for maps to apply. (See Chossat 
&; Iooss (1994) for a discussion of the centre manifold reduction in the similar problem of 
bifurcations from Taylor vortices.) Let ( be the eigenfunction associated with the critical FM, 
so that on the centre manifold, we can write 

A n = a„C + $K), (14) 

where 3> is the graph of the centre manifold. The unfolded dynamics takes the form 

a n+1 = g(a n ;fi), c n+1 = c n + (-l) n h(a n ; fi), (15) 

where g and h are the maps Q and h reduced to the centre manifold; g and h share the same 
symmetry properties as Q and h. 

In this paper, we only consider the case where T\ acts trivially. We therefore expect only 
generic bifurcations in the map g: saddle-node when FM = 1, period-doubling when FM = — 1 
and Hopf when there are a pair of complex FMs. The FMs for the full period map Mq, which 
are the squares of the FMs of g, will generically be either one or come in complex conjugate 
pairs. In particular, we do not expect Mq to have a FM = — 1; this mechanism for suppressing 
period-doubling bifurcations was discussed by Swift & Wiesenfeld (1984). 

Here we consider only the cases where g has a FM = +1 or —1. The normal form in the 
case FM = 1 is: 

a n+1 = + a n - a 2 n , c n+1 = c n + (-l) n h(a n ; fi), (16) 

to within a rescaling and a change of sign. The parameter [i is zero at the bifurcation point, 
and the fixed points of the a part of the map are a = when fj, is positive. The spatial 

translations are 

Co, ci = c + h(a; fi), c 2 = c , ... (17) 

We therefore have a c -parameterized family of solutions that vanish, in pairs, as fi is decreased 
through fx = 0. We interpret this bifurcation by considering the solutions with c = —\h{a\ fx), 
a = iy/Ti. In this case, we have a pair of pulsating wave solutions, translated with respect 
to the original PW by ^h(±^/fl; fi), which collide in a saddle-node bifurcation at fi = 0. The 
remainder of the family of solutions is obtained by translating this pair. 
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The case FM = — 1 is more interesting. The normal form in the supercritical case is 

a n +i = (-1 + v)a n ~ o? n , c„+i = c„ + (-l) n h(a n ;p), (18) 

with a fixed point a = and a period-two orbit a n = (—l) n ^/p. The dynamics of the spatial 
translations are 

c , ci = c + h(a Q ;p), c 2 = c + /i(a ; p) - h(-a ; p), 

(19) 

c 3 = c + 2h(a ;p) - h(-a ;p), ... 

Since /i(0;0) = 0, and generically ||(0;0) / 0, h(a ; p) and h(—a ;p) have opposite sign for 
small /i; this represents a symmetry-breaking bifurcation that leads to a solution that drifts 
along the group orbit of the PW. 

The main points of interest in this section are the approach that we have taken in analysing 
the instabilities of the group orbit of the spatio-temporally symmetric periodic orbit, and the 
observation that an instability of the pulsating wave with FM = 1 in the full-period map 
can lead to drifting solutions or not. Whether solutions drift can only be determined by 
examining the half-period map. In the next two sections, we apply our method to three- 
dimensional alternating pulsating waves and to alternating rolls, the latter having spatial as 
well as spatio-temporal symmetries. 



3. Three dimensions: alternating pulsating waves 

Alternating pulsating waves (APW) are the simplest three-dimensional analogue of the 
pulsating waves discussed in the previous section. These periodic oscillations have been 
observed in numerical simulations of three-dimensional compressible magnetoconvection with 
periodic boundary conditions in the two horizontal directions (Matthews et al 1995). They 
appear either after a series of global bifurcations (Rucklidge & Matthews 1995; Matthews et al 
1996) or in a Hopf bifurcation from convection in a square pattern (Rucklidge 1997), and are 
invariant under the combined operation of advancing one quarter period in time and rotating 
90° in space. 

The full symmetry group of the problem is the semi-direct product of the D 4 symmetry 
group of the square lattice and a two-torus T 2 of translations in the two horizontal directions, 
x and y. D 4 is generated by a reflection k x , a clockwise rotation by 90° p, 

k x : (x,y) -> (~x,y), p: (x, y) -> (y, -x), 
T a y. (x,y) — > (x + a (mod£),y + b (mod £)), 

where pr afi = n_ a p. 

As before, we assume that at p = 0, we have a known APW solution U 0i0 (x,y, z,t) with 
spatial period A in each direction and temporal period T; then f/ ,o satisfies: 

Uo.o(x,y, z, t) = pU (x, y,z,t + \T) = U Q (x, y,z,t + T) 

(21) 

= T Ai0 ^o,oOr, y, z, t) = T QtX U 0fi (x, y, z, t). 

There is a two-parameter continuous group orbit of APWs generated by translations: U a ^ = 
T a ,bUafi. We expand about this group orbit: 



U(x,y,z,t) = T Cx{t) , Cy{t) (U 0fi (x,y,z,t) + A(x,y, z,t)), 



(22) 
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where (c x , c y ) is a time-dependent translation around the group orbit and A is the perturbation 
orthogonal to the tangent plane to the group orbit. As before, we separate the evolution of 
the translations from that of the perturbation: 

d^4 dc dc 

— = g(A,U . ;p), -^ = h x (A,U 0fi ;p), -^ = h v (A,U , ;p), (23) 

where we keep track of the explicit time-dependence of Q, h x , and h y through the argument 
?7 0> o- The spatio-temporal symmetry of the basic state C/ ,o is manifest in Q, h x and h y as 
follows: 

G(pA, pU . ; p) = pG{A, U . ;p), 

h x (pA, pU 0fi ; p) = h y (A, U 0fi ; p), (24) 
h y (pA, pU 0fi ; p) = -h x (A,U afi ;p), 

where pU fi(t+ jT) = U , (t). It is convenient to introduce a complex translation c = c x + ic y 
and a corresponding h = h x + ih y , so pr c = T- lc p. 

As before, we define a time advance map acting on the perturbation so A(t) = Ai o (A(0)); 
this has the property 

Mlp = P M T T % + \ My = P *M T T % + \ Mlp 3 = p 3 M^/; 4 +t (25) 

because of the underlying spatio-temporal symmetry of the APW. The full period map Mq 
is then the fourth iterate of a map Q: 

Ml = P 4 M T 3T/4 Ml T/4 = P M T /4 p*Ml T/4 = (pM^y ee g\ (26) 

Instead of A4q , we consider Q = pAi^ 4 , which has no special properties under reflections and 
rotations, but does commute with T\, and r 1 \, which leave the underlying APW invariant. 

The dynamics of the perturbation is given by A n +i = G(A n ), where A(\T) = p 3 Ax, etc. 
Then 

,T/4 

Cl =co+/ h(Ml{Ao),U 0fi (t);p)dt = co + h(A ;p), (27) 
Jo 

where the map h = h x + ih y is invariant under the translations r A >0 and r A . After the next 
quarter period, we find 

r-T/2 



c 2 = ci + / h{M* T/i (A{\T)), U 0fl {t)-p) dt 

JT/4 
fT/2 

= ci+ / h(M t T/4 (p 3 Ai),U 0>0 (t);p)dt 

JT/i 

= Cl + / h{p 3 Mf T/ \Ai), p 3 U , (t - T/4); p) 
Jt/a 



T/4 
T/2 

. , , d/ 

T/4 

Cl +ih(Ai,p). (28) 



So the combined dynamics of the perturbation and the translation can be written as 

An+i = G(A n ;p), c„+i = c„ + i n h(A n ;p), (29) 
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where <?(0;0) = h(0;0) = 0. 

We consider the bifurcations of (29) only in the case where t\ i0 and t .\ act trivially. Note 
that, as in the case of pulsating waves, the generic bifurcations of APW are either steady state 
(FM = +1) or Hopf, since A4q = Q 4 . We consider bifurcations with FM = +1 of A4q only; 
generically, these occur when the linearization of Q has a FM of +1 or —1. Near a bifurcation 
point we reduce the dynamics onto the centre manifold 

a n +i = g(a„;n), c n+1 =c n + i n h(a n ; //). (30) 

When a FM = 1, once again we have a saddle- node bifurcation, this time involving pairs 
of APWs that are translated relative to each other. If a FM is —1, we have a n = (—l) n ^/pZ, 
and the spatial translations are: 

c , ci = c + h(a ;p), c 2 = c + h(a ; p) + ih(-a ; //), 

(31) 

c 3 = Co + ih(-a ;n), c 4 = c , ... 

This solution has no net drift (unlike in the two-dimensional problem), but travels back and 
forth different amounts in the two horizontal directions since, generically, h x (a ; p) ^ h y (a ; p). 
The solution remains invariant under advance of half its period in time combined with a 
rotation of 180°. To see this, we construct the solution U(x, y, z,t) at t = and t = \T using 
the solution in the c -parameterized family that satisfies c = — c 2 . Specifically, we insert the 
centre manifold solution A(0) = Aq = a a ( + $(a ), A(\T) = p 2 A 2 = p 2 (a Q ( + $(a )) in (22). 
We obtain 

17(0) =T Co ([/ ,o(0)+aoC + $(a )) 
U(±T) = r_ Co (£/ ,o(|T) + p 2 a ( + p 2 <S>(a Q )) 
= r_ co p 2 (C/ , (0)+a C + ^(a )) 

= P 2 U(0), (32) 

where we have suppressed the (x, y, z)-dependence of U, retaining only its t-dependence. 

Thus, in the simple case of APW, we cannot get drifting solutions in a bifurcation with 
FM = 1 for the time-T return map. We next consider the same bifurcation for the more 
complicated example of alternating rolls. This solution has the same spatio-temporal symmetry 
as APW but has extra spatial reflection symmetries. We shall see that in this case a particular 
symmetry-breaking bifurcation leads to two distinct types of drifting solutions. 



4. Additional spatial symmetries: alternating rolls 

Alternating rolls (AR) are created in a primary Hopf bifurcation from a k T 2 invariant 
trivial solution (Silber & Knobloch 1991). Like alternating pulsating waves, alternating rolls 
are invariant under the spatio-temporal symmetry of advancing one-quarter period in time and 
rotating 90° in space, but have the additional property of being invariant under reflections in 
two orthogonal vertical planes. Alternating rolls have been observed in three-dimensional 
incompressible and compressible magnetoconvection (Clune & Knobloch 1994; Matthews et al 
1995). 

For convenience in this section, we define p to be the combined advance of one quarter 
period in time followed by a 90° clockwise rotation about the line (x,y) = (0,0). Reflecting 
in the planes x = j\ or y = |A leaves alternating rolls unchanged at all times, so the sixteen- 
element group that leaves AR invariant is generated by k' x , k' v and p, where 

k' x : (x, y, z, t) -> (|A - x, y, z, t), 

n' y :(x,y,z,t) -> (x,±\ - y,z,t), (33) 
p: (x, y, z, t) -> (y, -x, z, t + \T). 
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Table 1. Summary of six types of bifurcations of alternating rolls, distinguished by the action 
of k x and n'y on the critical modes, and by the critical Floquet multipliers of Q. Isotropy 
subgroups (up to conjugacy) of the bifurcating solution branches are indicated, along with 
their order; in cases B+ and B— , there are two distinct solution branches. 



Case 


Action of k' x , K y on 
marginal modes 


Floquet 
multiplier(s) 


Bifurcation 
(drift or not) 


Isotropy 

subgroup (order) 


A i / i 1 \ 

A+(+l) 


ii/- /- 

•^x^yS — S> 


b M = + 1 


Saddle-node 
(no drift) 


(k x ,K y ,p) 


(IbJ 


A i ( 1 \ 

A+(-l) 


as A+(+l) 


t M = — 1 


Symmetry-breaking 
(no drift) 


v x ' y ' " / 


(8) 


A-f+1) 


k'k'„C = C 

K xC = K yC ~ ~C 


FM = +1 


Symmetry-breaking 
(no drift) 


(K x Ky , p) 


(8) 


A-(-l) 


as A-(+l) 


FM = -1 


Symmetry-breaking 
(no drift) 


(K x Ky, K X P~) 


(8) 


B+ 


C- = «iC+ = -«iC+ 


FM = ±1 


Symmetry-breaking 
(no net drift) 


(k' x ,P 2 ) 

(p) 


(4) 
(4) 


B- 


as B+ 


FM = ±i 


Symmetry-breaking 
(drift) 


Wy) 

Trivial 


(2) 
(1) 



The basic AR solution U .o(x,y, z,t) exists at (i = and satisfies 

U , (x, y, z, t) = pU 0fi (x, y,z,t+ \T) = U 0fi (x, y,z,t + T) 

= K' x U Qfi (x, y, z, t) = K' y Uo,o{x, y, z, t) (34) 
= T X , U 0fi (x, V, z, t) = T aA U 0fi (x, y, z, t). 

As in section 3, we expand about this basic solution and recover the map (29). The presence 
of extra reflection symmetries of the underlying solution manifests itself in the following way: 

G{k' x A) = k' v §(A), G{k' v A) = k' x G(A), 

hMA) = -h x (A), h x (K' y A) = h x (A), (35) 
h y (K' x A) = hy(A), h y {n' y A) = -h y {A). 

Note that the rotation in the definition of Q (26) implies that reflecting with n' x then applying 
Q is equivalent to applying Q then reflecting with n' y , since pK x = n' y p. In the terminology 

of Lamb & Quispel (1994), k' x and n' y are 2-symmetries of Q, that is, Q 2 (n' x A) = k' x Q 2 (A). 
In general, /c-symmetries arise when the spatial part of the spatio-temporal symmetry of a 
time-periodic solution does not commute with its purely spatial symmetries (Lamb 1997). 

The remainder of this section is devoted to the discussion of the codimension-one steady 
bifurcations of this problem. We do not consider bifurcations that break the spatial periodicity, 
so Ta,o and t ,a act trivially, nor do we consider Hopf bifurcations. The results are summarised 
in Table 1. 

We begin by noting that Q(K x K y A) = K,' x K y G(A), so n' x K y commutes with the linearisation C 

of Q, whereas k' x C = Cn' y . The eigenspaces of C are invariant under the reflection n' x K y . We 
assume the generic situation of one-dimensional eigenspaces, then each eigenfunction ( must 
be either even or odd under the reflection k' x k' v , i.e., k' x k' v Q = ( (case A) or k' x k' v ( = —( 

(case B), since (K x n y ) 2 is the identity. In case A, if ( = K x n' y ( is an eigenfunction of C with 
FM s, then k' x ( = k' v ( has the same FM: 



£«4C = K y £( = SK y ( = SK X (. 



(36) 
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Therefore, £ and n' x Q are linearly dependent; moreover, k' x 2 is the identity, so either k' x Q = C 
(case A+) or k' x ( = —( (case A—). Finally, these two cases are subdivided according to the 
value of the critical Floquet multiplier of C, (either +1 or —1) at the bifurcation point. 
Case B is rather different. Here we have k x C = —K y (, so 

£<C = KyCC = sk' v ( = -sk' x (. (37) 

Thus k' x C has FM = — s and is linearly independent of £, which has FM = s. We define £+ 
to be the eigenfunction of s and £_ to be the eigenfunction of —s, with = n' x C+ = — K ' y C+- 
There are two ways in which two Floquet multipliers s and — s can cross the unit circle: 
either at +1 and —1 (case B+) or at +i and — i (case B— ). Note that in the absence of the 
reflection symmetries these bifurcations would be codimension-two; here they occur as generic 
bifurcations. Since the FMs of the time-T map Mq are the fourth power of the FMs of Q, the 
effect of the symmetry in case B is to force a repeated FM = +1 in the map A4q . 
In case A, we write 

A n = a n ( + $(a„), (38) 

near the bifurcation point, where is the graph of the centre manifold. On the centre manifold 
we have A = k' x k' v A, so 

h x (A) = h x (K' x K' y A) = -h x (n' y A) = -h x (A) = 0, (39) 

where we have used (35). Thus in case A, h x and h y are identically zero, and no bifurcation 
will lead to drift along the group orbit of alternating rolls. 

The reflections k' x and K y act trivially in case A+. A FM = +1 leads to a saddle-node 
bifurcation of alternating rolls. The normal form in the case FM = — 1 gives a n = (— l)"a , 
from which the bifurcating solution U(t) can be reconstructed. Choosing the initial translation 
c to be zero, and suppressing the (x, y, z)-dependence of U, we have 

U(0) = U Ofi (0) + a ( + <S>(a ), U{\T) = p\U o , o (0) - a ( + $(-o )), 

U{\T) = p 2 (C/ o ,o(0) + a C + $(ao)), U{\T) = P {U^(Q) - a a ( + <D(-a )). ' 

Here it should be recalled that U 0t0 (^T) = p 3 U 0t0 (0), and that on the centre manifold 

A{\T) = p 3 A, = p 3 ( ai ( + $( 0l )) = p 3 (-a C + $(-ao)). (41) 

This solution satisfies 

U(t) = K ' x U(t) = K> y U(t) = p 2 U(t + ±T), (42) 

and thus has the same symmetries as "standing cross-rolls" , described by Silber & Knobloch 
(1991). 

In case A—, n' x and n' y act nontrivially, so the behaviour on the centre manifold is governed 
by a pitchfork normal form (a n = a ) when the FM = +1 and by a period-doubling normal 
form (a n = (— l) n a ) when the FM = —1. At leading order in a , the bifurcating solutions 
U (t) in the two cases are 

U(0) = *7 ,o(0) + aoC, U{\T) = p 3 (U , (0) ± a Q, 

U{\T) = p 2 ([/ , (0) + ooC), U{\T) = p(t/ ,o(0) ± aoC). 

These solutions are not invariant under k' x or K y (since these change the sign of Q, but are 
invariant under the product k x K y . In addition, U(t) = pU(t + \T) in the case FM = +1 and 
U(t) = K x pU(t + \T) in the case FM = -1. 
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Case B is more interesting. On the two-dimensional centre manifold, we write 

An = (-a n + b n )(+ + K + 6„)C- + *K, &„); (44) 
the form of this expression is chosen for later convenience. The map (29) reduces to 
(a n+ i,6„+i) = g(a n ,b n ;n), c n+1 = c n + i n (h x (a n , b n ; n) + ih y (a n , b n ; //)). (45) 
Since C- = <C+ = ~ K ' y C+, we have 

<A» = (a n + b n )( + + {-a n + b n )(- + <$(a„, &„), 

(46) 

re^At = (-a« - MC+ + ( a ™ ~ MC- + ^$(a„, &„); 

thus 

«4(a n , 6 n ) = (-a„, &„), «(,(a n , 6 n ) = (a n , -&„). (47) 
From this and from (35), we deduce that on the centre manifold 

hx^ni bn) hxi, ^ni b n ) ^i(^ni b n ) 

hy(& n , b n ) Hy{(lni fc n ) /iy( Q-n, &n)) 



(48) 



implying that h x (0, b; /x) = and /i a (a, 0; /i) = 0. Moreover, g inherits the symmetries (35) of 

K y g(a n , b n ) = g(K' x (a n , &„)), n' x g{a n , K) = g{K y {a n , b n )). (49) 
Thus the linearisation £ of g satisfies 

1 ^^f- 1 ?V (50) 



(51) 



-1 1 / ' 



which forces £ to be of the form 

a 

(3 



where a n , b n can be scaled so that a = 1. There is a bifurcation when /3 = +lor/3 = — 1, 
yielding FMs ±1 (case B+) or ±i (case B— ), respectively. 

In order to analyse the dynamics near the bifurcation point, we compute the normal form 
of the bifurcation problems, expanding g as a Taylor series in a and b. The reflection symmetry 
k' x n'y prohibits quadratic terms, and all but two of the cubic terms can be removed by near- 
identity transformations. We thus have the unfolded normal form, truncated at cubic order, 
in the two cases B+ and B— : 

O-n+1 = b n , 

b n +i = ±(1 + fi)a n + Pal + Qa n b 2 n , (52) 
c n +i = c n + \ n {h x (a n ,b n ;n) + ih y (a n , b n ; //)), 

where (i = at the bifurcation point and P and Q are constants. 

Lamb (1996) deduced the (a, b) part of this normal form for a local bifurcation of a map 
with a Z 2 x Z 2 2-symmetry group, appropriate to the case under study here. We have chosen 
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our scalings and near-identity transformations to match Lamb's notation. Lamb described the 
period-one, two and four orbits that are created in the bifurcation at p = and calculated 
their stability as a function of the constants P and Q; we will interpret those results in terms 
of bifurcations from, and drift along, the group orbit of alternating rolls. 

In case B+, there are three types of orbits created. The first is a period-two orbit 
(a ,0) <-> (0, a ), with = p + Pa§. From this and the symmetries (48) of h, we deduce 
the drift of the solution at each iterate: 

Co, ci = c + h x (a Q ,0;p), c 2 = c + h x (a , 0; p) - h y (0, a ; p), 

( 53 ) 

c 3 = c - /j y (0, a ; ju), c 4 = c , ... 

There is no net drift along the group orbit in this case. Moreover, the c -parameterized 
family of solutions drift to and fro in the x direction only since c„ — c is real. Consider 
c = ^(—h x (a ,0; p)+h y (0, a ; p)) = — c 2 , where c is real and thus corresponds to a translation 
in the x-direction. The reconstructed solution U(t), at leading order in a , satisfies 

17(0) = T Co ([/ ,o(0)-aoC+ + aoC-), U{\T) = T Cl p 3 ([/ . (0) + a C+ + a„C-), , . 

(54) 

17(§T) = r_ co p 2 ([/ ,o(0) - a C+ + a C-), I7(f T) = r_ Cl/ 9([/ , (0) + a C+ + a C-), 

so we have U(t) = K y U (t) = p 2 U(t + \T). The conjugate orbit, (0, a ) <-> (a , 0), has symmetry 
(k' x ,p 2 ) and does not drift at all in the x direction. 

The second and third types of orbit created in case B+ are a period one orbit (a , a ) and 
a period two orbit (a , — a ) <-> (— a , a ), with = /i + (P + Q)ag i n both cases. These orbits 
are mapped to each other by k' x or by K y , so we consider only the orbit with period one. The 
translations at each iterate are 

Co, ci = c + h(a ,a Q ;p), c 2 = c + h(a , a ; p) + ih(a , a a ; fj,), 

( 55 ) 

C3 = c + i/i(a ,a ;/i), c 4 = c , ... 

This orbit also has no net drift, and by choosing c = — 1(1 + i)^(«o, Oo; A*), w e have Ci = ic . 
The reconstructed solution U(t), at leading order in a , satisfies 

U(0) = r Co ([/ ,o(0) + 2a C-), ?7(iT) = T ico/ 9 3 (?7o,o(0) + 2a C-), / x 

(56) 

U{\T) = t_ CoP 2 ([/ ,o(0) + 2a C_), U{\T) = r_ ico/ 9([/ , (0) + 2a C-), 

so C7(t) = pU(t+jT), and the isotropy subgroup is (p). This solution has the same symmetries 
as the alternating pulsating waves described in section 3, so alternating pulsating waves may 
be created in a symmetry-breaking bifurcation of alternating rolls. The period-two orbit 
(a ,— a ) <-> (— a ,a ) has the conjugate isotropy subgroup (k' x k' v p). 

Finally, we turn to case B— . Here, there are two types of periodic orbit created in 
the bifurcation at p, = 0, and in this case they are both of period four. The first orbit is 
(a , 0) — > (0, — a ) — > (— a , 0) — ► (0, a ), with = — p + Pal in (52). The translations are 

Co, ci = c + h x (a a ,0;p), c 2 = c + h x (a , 0; /x) + h y (0, a ; p), 

(57) 

c 3 = c + 2h x (a Q , 0;p) + h y (0,a ;p), c 4 = c + 2h x (a , 0; //) + 2h y (0,a Q ; p), . . . 

Note that c„ — c is real so there is no drift at all in the y direction, but there is a systematic 
drift in the x direction. The reconstructed solution U (t) satisfies 

17(0) =T Co ([/ , (0)-aoC+ + aoC-), U{\T) = r Cl p 3 (f/ ,o(0) - a C+ - o C-), 

( ) 

f/(iT) = r C2 p 2 ([/ ,o(0) + a C + - a C-), U{\T) = r C3 p(U 0fi (0) + a ( + + a C-), 
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so we have U(t) = n' y U(t). A conjugate orbit, started a quarter period later, has isotropy 
subgroup (k' x ) and drifts systematically in the y direction. 

The second type of orbit created in case B— is (do, do) — > (a ,— a ) — > (— a ,— a ) — > 
(— a , a ), with = — /i + (P + Q)oq. The translations are 

c , ci = c + h x (a a ,a ;p) + ih y (a ,a ; p), 

c 2 = c + (1 + i)(/i x (a ,a ;/i) + h y (a ,a ; p)), 

(59) 

C3 = c + (2 + i)h x (a , oo;/i) + (1 + 2i)/ij, (a ,a ; //), 

c 4 = c + (2 + 2i)(/i x (a,a;/x) + /ij,(a , a ; //)), ... 

This corresponds to a solution that drifts along the diagonal, with a wobble from side to side 
as it goes. At leading order in a , the reconstructed solution U(t) satisfies 

Lf(0)=T Co ([/ ,o(0)+2a C-), U{\T) = r Cl p 3 ([/ ,o(0) - 2a C+), , N 

(60) 

U{\T) = r C2 p 2 ([/ ,o(0) - 2a C-), U{\T) = T C3 p{U 0fi {0) + 2a C+), 

which has fully broken the spatial and spatio-temporal symmetries of the underlying 
alternating rolls solution. 

We close our discussion of case B— by considering spatio-temporal symmetries of the 
drifting solutions that are obtained by moving to an appropriate travelling frame. These 
symmetries are not listed in Table 1. We first consider the solution (58) that drifts in the y 
direction; it has the following spatio-temporal symmetry 

U{t) = K' xP 2 T Ca _ C2 U(t + \T). (61) 

Next we consider the solution (60) that drifts along the diagonal; it has spatio-temporal 
symmetry 

U(t)=K' y pr ic ,. Cl U(t+\T). (62) 

In summary, we have examined the six different cases in which alternating rolls undergo a 
bifurcation with FM = +1 in the full period map. All six bifurcations preserve the underlying 
spatial periodicity of the alternating rolls, but may break the spatial and spatio-temporal 
symmetries. The 2-symmetry present in the B cases forces two Floquet multipliers to cross the 
unit circle together, and we find two branches of bifurcating solutions, with distinct symmetry 
properties. It is only in case B— , with Floquet multipliers ±i in the map Q, that the bifurcation 
leads to systematically drifting solutions: one solution drifts along a coordinate axis, while the 
other drifts along a diagonal. 

We finish this section by presenting an example of a bifurcation from alternating rolls to 
a drifting pattern, which we interpret as an instance of a B— bifurcation. We have solved 
the PDEs for three-dimensional compressible magnetoconvection in a periodic 2x2x1 box, 
using the code of Matthews et al (1995). The PDEs and description of the parameters and 
numerical method can be found in that paper. Figure 2 shows an example of an alternating roll 
at times approximately 0, jT, |T and |T; the two reflection symmetries k' x and n y in planes 
x = j\ and y = jX and the spatio-temporal symmetry of advancing a quarter period in time 
followed by a 90° rotation about the centre of the box are manifest. Increasing the controlling 
parameter, the temperature difference across the layer, leads to the solution in Figure 3: the 
data are shown at times 0, jT, |T, |T, T and 2T. The only spatial symmetry remaining is 
the invariance under n' , and the spatio-temporal symmetry has been broken. By comparing 
frames (a), (e) and (f) at times 0, T and 2T, it can be seen that the solution is drifting slowly 
leftwards along the x-axis. Moveover, a drift symmetry n' x p 2 T Co _ C2 conjugate to (61) can be 
seen by comparing frames (a) and (c). All evidence points to the bifurcation being of type 
B— (though we have not computed the Floquet multipliers and critical eigenfunctions in the 
PDEs). 
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(a) (b) 




x X x X 

Figure 2. Alternating rolls in three-dimensional compressible magnetoconvection, starting 
with parameter values from Matthews et al (1995). The four frames are (approximately) 
at times (a) t = 0, (b) t = \T, (c) t = \T and (d) t = §T. The frames show contours 
of the vertical velocity in a horizontal plane in the middle of the layer: solid lines denote 
fluid travelling upwards, dashed lines denote fluid travelling downwards, and the dotted line 
denotes zero vertical velocity. The spatial symmetries k' x and k' v are manifest, as is the spatio- 
temporal symmetry of advancing a quarter period in time followed by a 90° rotation (counter- 
clockwise in this example). The dimensionless parameters are: the mid-layer Rayleigh number 
(proportional to the temperature difference across the layer) R = 2324; the Chandrasekhar 
number (proportional to the square of the imposed magnetic field) Q — 1033; the Prandtl 
number a = 0.1; the mid-layer magnetic diffusivity ratio £ = 0.1; the adiabatic exponent 
7 = 5/3; the polytropic index m = 1/4; the thermal stratification 9 = 6; the mid-layer plasma 
beta [3 = 32; and the horizontal wavelengths A = 2 in units of the layer depth). 



5. Conclusion 

We have developed a technique for investigating the possible instabilities from continuous 
group orbits of spatio-temporally symmetric time-periodic solutions of partial differential 
equations in periodic domains. Our approach is based on centre manifold reduction and 
symmetry arguments. It is in the spirit of earlier work by Iooss (1986) on bifurcations 
from continuous group orbits of spatially symmetric steady solutions of partial differential 
equations. We have treated three examples that arise in convection problems: pulsating waves 
in two dimensions, and alternating pulsating waves and alternating rolls in three dimensions. 




Figure 3. After a bifurcation of type B— , the alternating rolls begin to drift. The parameter 
values are as in Figure 2, but with a higher thermal forcing: R = 3000 and Q — 1333. The 
frames are (approximately) at times (a) t = 0, (b) t = \T, (c) t = \T, (d) t = §T, (e) t = T 
and (f) t — 2T. Note how all spatial and spatio-temporal symmetries have been broken, with 
the exception of k v , a reflection in the plane y — \\ (modulo a slight shift in the periodic box). 
The slow leftward drift of the pattern can be seen by comparing frames (a), (e) and (f), In 
addition, a drift symmetry K' x p 2 T C0 - C2 , conjugate to (61), can be seen by comparing frames (a) 
and (c) or (b) and (d). 
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A simple bifurcation can lead to drifting solutions in the case of pulsating waves but not 
alternating pulsating waves. The additional spatial symmetries of alternating rolls can force 
two Floquet multipliers to cross the unit circle together; this degeneracy can lead to drifting 
solutions, as in the numerical example presented in the previous section. We have related our 
work to the theory of fc-symmetries developed by Lamb & Quispel (1994). 

Our approach can readily be applied to other problems. In the future, we plan to tackle 
spatial period doubling and multiplying, where the r A symmetries do not act trivially; such 
instabilities are relevant to simulations of convection carried out in larger boxes (Weiss et al 
1996), and will be related to the study of the long- wavelength instabilities of alternating rolls 
(Hoyle 1994). We also plan to examine the case of the hexagonal lattice: a Hopf bifurcation 
on a hexagonal lattice leads to a wide variety of periodic orbits with different spatio-temporal 
symmetries (Roberts et al 1986). Finally, we plan to investigate the effect of including the 
extra Z 2 mid-layer reflection symmetry that arises when making the Boussinesq approximation 
for incompressible fluids. 
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